Mon. Not. R. Astron. Soc. 000, (2005) Printed 2 February 2008 (MN JATfcjX style file v2.2) 



The core of the Canis Major galaxy as traced by Red 
Clump stars 

M. Bellazzini 1 *, R. Ibata 2 , N. Martin 2 , G.F. Lewis 3 , B. Conn 3 , M.J. Irwin 4 

l^"} i 1 INAF - Osservatorio Astronomico di Bologna, via Ranzani 1, 40127, Bologna, Italy 

f"*") , 2 Observatoire de Strasbourg, 11, rue de I'Universite, F-67000, Strasbourg, France 

f^*^ ■ 3 Institute of Astronomy, School of Physics, A29, University of Sydney, NSW 2006, Australia 

' 4 Institute of Astronomy, Madingley Road, Cambridge, CB3 0HA, U.K. 

> ' 

o 

Accepted for publication by MNRAS, 30 November 2005 

o 

CO 

ABSTRACT 

CNI ■ The recently-discovered stellar system in Canis Major is analysed using He-burning 

Red Clump stars as tracers. Canis Major turns out to be the strongest and most 
spatially confined overdensity of the whole Galactic Disc, both in terms of number 
density and of statistical significance. When projected onto the Galactic Plane, it 
appears as an elongated and compact overdensity extending from / ~ 200° to I ~ 280° 
with a roundish core toward I ~ 240°. We find that the main body of the system 
' has an integrated absolute magnitude My = —14.4 ± 0.8, a central surface brightness 

f*^ , fiy.o — 24.0 ±0.6 mag/arcsec 2 and a line-of-sight profile peaked at Dq = 7.2 ± 1.0 kpc 

with Half Width at Half Maximum ~ 2.0 kpc, in excellent agreement with the results 
obtained with widely different tracers (M giants and Main Sequence stars) in previous 
analyses. The mean distance to the main body of Canis Major is observed to increase 
with increasing Galactic longitude, from Dq ~ 6.3 kpc at I ~ 225°, to Dq — 9.3 kpc 
Jjj | at I ~ 265°, in good agreement with the predictions of our more recent N-body 

simulation that models CMa as a dwarf galaxy being accreted in a planar orbit onto 
the disc of the Milky Way. We confirm that the Canis Major system has all the 
characteristics of the relic of a dwarf galaxy seen on top of a large-scale overdensity 
that we detect all over the third and fourth Galactic quadrants (180° ^ I ^ 360°, 
with a strong maximum around I — 290° and b > —5°) that is identified as the stellar 
^ , component of the southern Galactic Warp. On the other hand, the possibility that a 

peculiar deformation/asymmetry of the outer Galactic Disc may be at the origin of the 
observed distribution of overdensities toward CMa cannot be definitely ruled out with 
the data presented in this paper. We also address a recent claim that Canis Major is 
on the outskirts of a larger "Argo" structure centred at I ~ 290° . Our analysis shows 
that the stellar populations in the latter are distributed over a very large distance 
range along the line of sight, and do not give rise to a narrow peak in density in this 
direction, contrary to what is observed in Canis Major. This suggests that the Argo 
structure is likely due to Galactic asymmetries such as the Warp. 
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1 INTRODUCTION 

During the last decade, it has been fully realised that the 
relics of the latest accretion events that contributed to the 
assembly of the Milky Way can be observati onally iden- 
tified and studied in deta i l (see, for example [ ibata et alJ 
1200 ll: iNewberg etafl 120021: iMaiewski et alJ l200l and ref- 
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erences therein). These findings simultaneously provide a 
qualitative confirmation and a new challenge for the cur- 
rent cosmological models in which the growth of large 
galaxies is d riven by the process of hierarchical merge r 
of sub-units JWhite fc Reeslll97& IWhite fc Frenkl Il99lll . 
The tidal disruption of dwarf galaxies within the Galactic 
potential may lead to the production of long-lived stellar 
streams (as in the ca se of the dwarf spheroidal Sagittarius 
galaxy - Sgr dSph, see llbata fc Lewisll998l : llbata et al.l200ll : 
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Ibata et al. 1 12001 llvezic et all Eoocl; iNewberg et al l EoO 
Maiewski et alJ 120031: iBellazzini. Ferraro fc Ibatal 12003* 

that may reveal fundamental information about the pro- 
cess of disruption, the mass distribution within the Galac- 
tic halo of Dark Matter, its degree of dumpiness, etc. (see 
Ibata et al.ll200ll: iHelmil 120041: I Johnston. Law fc Maiewski 



2003: lLaw. Johnston fc Maiewskill2005h . 

The recent discovery of a large stellar stream 
nearly coplanar with the Galactic disc (the Canis Ma- 
or/Monoceros Ring ; hereafter the Ring, for brevity 



jor/Monoccros rung ; hereafter the King, tor b revity 

Newberg et al.l|2002t"[Yannv et al."1l2003t llbata et alJl2003t 
2003tlRocha-Pinto et alll2003l : ICrane et al] 



Maiewski et al. 

20031: IConn et al.1 12005a 1 la) suggests that the accretion of 
small stellar systems may have had a considerable role also 
in the assembly of the disc components of the Milky Way 
(in particular, the Thick Disc) as envisaged by the latest de- 
tailed m odels of disc-galaxy fo rmation within a c osmological 
contex t llAbadi et alJl2003allbl: iHelmi et al.ll2003h . 

In iMartin et alJ (l2004al . hereafter Pap-I) we reported 
on the identification of a possible new stellar relic (the 
Canis Major dwarf galaxy, hereafter CMa) located at 7 — 
8kpc from the Sun whose approximate centre lies around 
(Z; 6) ~ (240°; —7°), and that may be the progenitor of th e 
Ring (see also iPenarrubia et al-lfiooit iDinescu et al.l [2005). 
CMa was identified as a strong elliptical-shaped overden- 
sity of M-giants by the comparison of star counts in North- 
ern and Southe rn Galactic hemis pheres (from the 2MASS 
All Sky Survey. ICutri et al.ll2003T) . The structure was sug- 
gested to be possibly a s sociated with some globular (bu t see 
IPenarrubia et ajj |2j)05jjMartincz-Dclgado_ej L _aljj2j^ 
open clusters ijMartin et al j l2004at IBellazzini et alJ l2004a ; 
ICrane et al] l2003h . The relic appears to have luminosity 
and size similar to the Sgr dSph (Pap-I). T he optical Color 
Magnit ude Diagram (CMD) we obtained in IBellazzini et al] 
(2004a, hereafter Pap-II) revealed a narrow and well-defined 
Main Sequence typical of an intermediate to old (age~ 4—10 
Gyr) and moderately metal deficient ([M/H] ~ —0.5) stel- 
lar system; these results have been fully confirmed by th e 
deeper CMD presented bv lMartinez-Delgado et all ]2005al) . 
A Blue Plume of possibly young stars or blue stragglers - also 
typical of dwarf sp heroidal galaxies - ha s been detected in 
both optical CMDs. ICarraro et al] J2005I) detected the Blue 
Plume population in the background of several open clusters 
in the third Galactic quadrant. While these authors conclude 
that this population should be associated to the Norma and 
Perseus spiral arms of the Milky Way, it turns out that all 
the considered fields lie within the region of the sky where 
CMa is detected (from their Table 1 219° sC I < 254°; 
—6.3° ^ b ^ +1.8°, see Pap-I and Sect. 3, below), and in the 
same range of distances as CMa (6.0 kpc ^ Dq ^ 11.7 kpc, 
see Pap-I, Pap-II and Sect. 4, below). Moreover, t he Blue 
Plume populations studied bv lCarraro et all ll2005l) appear 
to follow a distance - Galactic longitude trend quite similar 
to that found here for the Red Clump population of CMa 
(see Sect. 4.1.1, below). 

The first results of a large spectroscopic survey pro- 
vided an estimate of the systemic radial velocity and ve- 
locity disp_M|sjoiiof_CMa (V r — 110 km/s; a — 13.0 
km/s iMartin et all l2004bL hereafter Pap-III), while an es- 
timate of the systemi c proper motion has b een obtained by 
iMomanv et alj (120041) . IMartin et al.1 feOOSl hereafter Pap- 



Branch and Red Clump stars around the centre of CMa, 
confirming the essence of the results of Pap-III, namely 
that CMa stars have a systemic velocity and a velocity dis- 
persion quite different from what expected from Galactic 
disc stars at the considered distance and quite typical of 
Galactic dwarf satellite galaxies. Moreover they detected a 
clear distance - radial velocity gradient among CMa stars 
that can be explained as the effect of on-going tidal dis- 
ruption of the stellar system. A significant refinement of 
the pro per motion est i mate was also obtained in Pap-IV. 
Finally. iDinescu et al] l)2005S) very recently obtained a new 
accurate estimate of the proper motion of CMa and con- 
cluded that (a) the system has a significant component of 
motion perpendicular to the Galactic Plane (W ~ —50 
km/s) and that (b) this W component of the space mo- 
tion of CMa is in the opposite direction with respect to 
the expected W motion of stars associated with the Galac- 
tic Warp. Other interesting investigations of the CMa relic, 
not strictly related to the present work or in a too prelimi- 



nary stage to dc iuny usciui ior t 
in iKinman. Saha fc Pied J2004D: 


nc analysis, can dc iouna 
Forbes. Stradcr & Brodid 


(2004j):ISbordone et al] 120051) : iMateu et alJ 


(2005). 



IV) analysed the kinematics of more than 1400 Red Giant 



density may be an effect of the stel l ar component of the 
Galactic Warp JPiorgoyski fc Sosinl Il989l lYusifovl 12004 
lLopez-Corredoira et al] l2002t IVig. Ghosh fc Qihal 120051) . 

whose broad southern maximum (as a South/North den- 
sity asymmetry) lies in the range 240° ^ I ^ 3 00° 
jLopez-Corredoira et all20ollVig. Ghosh fc Oihal2005l) . In 
Pap-II and Pap-III, in particular, we showed that such a 
conclusion can be reached only if the distance distribution 
of the adopted tracers (M giants) is neglected and that CMa 
is in fact an additional overdensity of different nature with 
respect to the large scale Warp of the Galactic disc. On the 
other hand the disentanglement of the two structures may 
prove difficult, limiting our possibility to obtain a detailed 
global description of CMa and the Warp (see Pap-III, for 
this discussion) . A global analysis of CMa is also hampered 
by the large degree of contamination from Galactic stars, by 
the high and variable extinction toward the Galactic plane 
and by its proximity to us that imply a very low surface den- 
sity in star counts per unit solid-angle of sky (see Sect. 4.4, 
below) . 

It is clear that a large scale study of CMa using a differ- 
ent tracer than the M-giants used in the previous wide-field 
studies of the system would provide a very useful indepen- 
dent check of the results of Pap-I and may help us to un- 
derstand the relation between CMa and the Warp. In Pap-I 
we showed that the Red Clump (RC) of Helium-burning 
stars of the CMa system is clearly detected in the 2MASS 
database and in Pap-II we used the RC stars to show that 
CMa is an independent structure, superposed on the Galac- 
tic Warp. Here we extend the use of RC stars to trace the 
whole structure of the outer Galactic Disc, focusing on the 
CMa relic and of the Galactic Warp in its surroundings. In 
Sect. 2 we describe the adopted dataset and assumptions, 
and we introduce our analysis method. In Sect. 3 we study 
the spatial distribution of overdensities rising from asym- 
metries between the Southern and Northern Galactic hemi- 
sphere, focusing on the relation between Canis Major and 
the Galactic Warp. In Sect. 4 we derive some fundamen- 
tal parameters of CMa as a stellar system, and, finally, in 
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Figure 1. Left panel: Infrared Colour Magnitude Diagram of 
a 6° X 5° field near the centre of the CMa overdensity. The Red 
Clump of CMa is indicated by the arrow. The vertical lines enclose 
the colour range in which we select candidate RC stars. Right 
Panel: Luminosity Function of the selected RC stars (dotted black 
line). Note the strong bump associated with the RC of CMa at 
K ~ 13.0. The Luminosity Function obtained exactely in the 
same way from the corresponding synthetic sample extracted from 
the R03 model (see Appendix A) is plotted as a red continuous 
line. The synthetic LF has been rescaled by a factor 1/3 to match 
the total number od stars in the two samples. 

Sect. 5 we summarise and discuss the main results of these 
analyses. 



2 RED CLUMP STARS AS TRACERS OF 
OUTER DISC STRUCTURES. 

We extracted from the 2MASS Point Source Catalogue 
the sources comprised between I = 40° and / = 320°, 
having —30° ^ b ^ 30°. The selected range of longi- 
tude and latitude cover the largest part of the Galactic 
Disc, avoiding only the surroundings of the Galactic Cen- 
tre, where the stellar crowding and the interstellar ex- 
tinction may prove too severe for a useful analysis in the 
present context. Stars were selected according to the qual- 
ity of the photom etry and to avoid spurio us and contami- 
nated sources as in Bellazzini et al 1 J2003bl) . Moreover, only 
stars having photometric uncertainties < 0.3 mag both 
in J and K magnitude have been included in the sam- 
ple. The magnitudes of all the stars were corrected for 
extinction using the c olour excesses interpolated on the 
COBE/DIRBE maps dSchlegel. Finkbeiner fc David ll99Sl) 
and modifie d according totot^ ^^^o ^Mana^^Beeial^O^I . 
as done in iMomanv et all il2004Tr and Pap-III. Obviously 
a 2-dimensional extinction map cannot account for the 3- 
dimensional distribution of the intersellar dust, thus the 
adopted correction may not always be good for nearby 
stars, but should be quite appropriate for the distant Red 



Clump stars we a r e inte rested in. The extinction laws by 
iRieke fc Lebofskvl lll985tl have been adopted. In the fol- 
lowing J, H, and K magnitudes denote reddening-corrected 
magnitudes Jo, Ho, and Ko, unless otherwise stated. Only 
stars with K < 15.0 were retained for the analysis. 

In Figure 1 we show the J-K, K CMD of a low-extinction 
field near the centre of the CMa structure. The diagram is 
dominated by Main Sequence (MS) stars of the Galactic 
disc for J — K < 0.5, for J — K > 0.8 the M giants are con- 
tributed by both the disc and CMa, and the intermediate 
colour strip, 0.45 ^ J — K ^ 0.70, enclosed by continuous 
lines in Fig. 1, should be dominated by RC stars located at 
various distances. While the M giant sequence of CMa can 
be revealed only with the subtraction of a control field sym- 
metric with respect to the Galactic Plane (Pap-I) the Red 
Clump of the system is clearly visible without subtraction 
in the range 12 ^ K ^ 14 with an apparent peak around 
K — 13, as readily shown by the histogram in the left panel 
of Fig. 1 (dotted black line), in excellent agreement with 
Pap-I. Hence RC stars are very promising tracers of the 
CMa structure. In the following discussion we will limit our 
analysis to the stars in the colour strip enclosing the typical 
RC stars, 0.45 ^ J — K 0.75, as in Pap-II (e.g., the region 
enclosed within parallel continuous lines in Fig. 1). In the 
left panel of Fig. 1 we also plotted the Luminosity Function 
(LF) obtained with the same selection criteria from a syn- 
thetic sample of the same region o f the sky extracted from 
the Galaxy model bv lRobin et alJ <l2003l . hereafter R03, red 
continuous line; see Appendix A for details). The R03 model 
includes a warped and flared disc (whose parameters have 
been fixed according to observations) and the effect of the 
displacement of the Sun with respect to the Galactic Plane 
by 15 -16 pc (toward the Galactic North. iHammerslev et all 
1995, R03). It is interesting to note that the R03 model 
reproduces excellently the observed LF for K < 12.3 and 
K > 13.7, but completely lacks the strong bump around 
K — 13 that we attribute to the RC of Canis Major. 

With respect to M giants, RC stars may suffer from a 
larger degree of contamination by unrelated sources. While 
for J — K ^ 1.0 the only possible contaminants for M gi- 
ants are M dwarfs at faint magnitudes (K 12.0), our 
colour-selected RC sample may be contaminated by K gi- 
ants essentially at any magnitude, and significant num- 
bers of both M and G dwarfs can be pushed into the 
RC colour range by photome tric errors for K > 14.0 (see 
lLopez-Corredoira et alj|2002l and Appendix A, for discus- 
sion) . However all the results presented below are in the form 
of subtractions of star counts in the Southern and North- 
ern Galactic hemispheres. Since the degree of contamination 
from Galactic sources as a function of magnitude should be 
similar in both hemispheres, the subtraction should min- 
imise any spurious effect associated with contaminants (see 
Appendix A, below). Moreover, all the results presented be- 
low remain unaltered if we limit our sample to RC stars 
having K ^ 14 .0, where the contamination b y dwarfs is not 
overwhelming tLopez-Corredoira et alJl2002l) . In Appendix 
A we study the problem of contamination by dwarfs in some 
detail, showing that (a) in spite of the large fraction of dwarf 
that may contaminate our samples and of the presence of 
North-South differences in the degree of contamination, the 
adopted technique traces trustworthy real asymmetries in 
the density of distant giants, for K ^ 14.0, and (b) the in- 
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Figure 2. Map of the extinction (in K band) for the 
range of Galactic longitude and latitude considered in 
the present analysis. The COB E-DIRBE reddening maps 
ISchlegel. Finkbeiner fc David 11998 ) have been interpolated on 
a grid whose knots are spaced by 0.5 deg both in I and b. 



trinsic asymmetry due to the ~ 15 pc displacement o f the 
Sun above the Galactic plane llHammerslev et "al]ll995fl has 
a negligible effect on our analysis. 

On the other hand RC stars should be intrinsically more 
numerous than M giants and they provide a completely 
independent and more reliable distance sc ale compared to 
the photometric parallax of M giants (see iBellazzini et alJ 
j2003bl) . lMaiewski et ail <2003ft . Pap-I, Pap-II and Pap-III). 

If not otherwise stated, we will use in the following 
analysis only stars with |6| > 5.0° to avoid the regions 
where the extinction is most severe. Fig. 2 shows that once 
the |6| > 5.0° strip is excluded, the region considered in 
the present analysis is largely free from high extinction re- 
gions (e.g. where Ak > 0.3 mag). Since most of the struc- 
tures we are interested in are enclosed in the latitude strip 
|6| ^ 15.0°, most of the following analysis is limited to this 
range. Note that the 2M ASS Point Source C atalogue is com- 
plete down to K$ = 14. 3 jCutri et all2003fl . hence the lumi- 
nosity range in which most of the RC stars we will consider 
lie ( 11.0 < K < 14.0) is not expected to suffer from prob- 
lems of incompleteness. 



2.1 CMa and the Warp as South-North 
overdensities 

It has to be noted that the origin of disc warps and the rela- 
tionshi p between gas e ous and stellar warps is far from under- 
stood fainnevllT99l iKuiiken fc Garcia-Ruizl 1200 ll) . More- 
over the characterisation of the stellar Warp of the Milky 
Way is quite poor and uncertain, mainly because of the 
problems inherent with the analysis of stellar fields at low 
Galactic latitudes. Finally, it is important to realise that 
the adopted parametrisation may describe different char- 



acteristics of the warp and ma y induce some confu s ion in 
the use of terms. For instance, iDiorgovski fc Sosinl (ll989T> 
identified the stellar warp as a sinusoidal variation of the 
mean latitude of their tracers (namely IRAS sources) as a 
function of longitude. In this case the "maximum of the 
warp" coincides with the maximum angular height above 
the Galactic equator, similar to the parametrisation adopted 
when dealing with the gaseous warp. On the other hand, 
lLopez-Corredoira et alJ i2002l) described the warp as the 
sinusoidal variation of the ratio of star counts in North- 
ern and Southern Galactic hemisphere as a function of 
longitude, hence their "maximum of the warp" is the re- 
gion where the Southern (or Northern) projected overden- 
sity reaches its maximum. The two definitions may pro- 
vide very different views, depending on the structure of 
the disc and on the relative 3-D positions of the observed 
structures and the Sun. Since the CMa galaxy was iden- 
tified as a South/North overdensity it should be com- 
pared with the Warp as a S outh/North overdensity (a la 
lLopez-Corredoira et aljl2002l) , since this is the only relevant 
parameter in this context. The parametrisation of the am- 
plitude of the Warp as the artificial shift to apply to the 
Galactic plane to match the star co unts in the Northern 
and Southern hemisphere (as done bv lMomanv et al"ll2004h 
is not a fair description of South/North overdensities in the 
present context, since it turns into a comparison of stars ly- 
ing at different distances, as shown in Pap-III. Hence, in the 
following we will always deal with CMa and the Warp as 
South/North overdensities. 

2.2 The construction of subtracted density maps 

Most of the following analysis will make use of subtracted 
density maps obtained form the above described sample of 
color-selected RC stars. In this section we will describe in 
detail how these maps are obtained. 



example [Paczvnski & Stanek 


Il998l: ISaraiedini et al 


2002 


Babusiaux & Gilmorel 12005: 


Lopcz-Corrcdoira et al. 


2002 



and references therein). In Pap-II we assumed Mk = — 1.5± 
0.2 for the RC of CMa, since , acco rding to the theoretical 
models by ISalaris fc GirardU i2002ft . this interval encloses 
the absolute K magnitude of the RC for populations of 
metallicity -0.4 < [M/H] < -0.7 and age between 4 Gyr 
and 10 Gyr, i.e. the proper age and metallicity for CMa 
(see Pap-II). Here we adopt the same assumption to obtain 
the distance of each RC star in our sample. These individ- 
ual distances should be regarded with some caution since 
the adopted distance scale is optimised for the RC of CMa 
and may provide misleading answers when is applied (a) 
to contaminant sources that are not RC stars (see sect. 2, 
above), and (b) to RC stars outside the age and metallic- 
ity ranges outlined above. On the other hand, the range 
of age and metallicity of stellar populations in the outer 
Galactic Disc that produce RC stars should not be too dif- 
ferent from that outlined above. Hence, the assumed global 
distance scale is expected to provide a reasonable descrip- 
tion of the overall spatial distribution of RC stars in this 
region of the Galaxy and a reliable representation of the 
CMa system. In particular, it should be considered that RC 
stars are much less sensitive to metallicity and age varia- 
tions in the considered populations compared to M-giants 
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(see. IMaiewski et al.ll2003l: ICrane et alJl2003t IMartin et"al] 
l2004al : iRocha^mtoetall 2005L and references therein) . Fi- 
nally, it should also be recalled that in the following we will 
deal only with subtracted density distributions: the subtrac- 
tion should cancel out many of the above effects if, as it very 
reasonable to assume, at a given Galactocentric distance the 
age and metallicity distributions are similar in both Galactic 
hemispheres. 

We used the distances of individual RC stars in our 
sample to obtain their Cartesian Heliocentric coordinates: 



x — DQCOs(Z)cos(6)[kpc] 
y = DQsin(l)cos(b) [kpc] 
z = Z?Qsm(&)[kpc] 



(1) 
(2) 
(3) 



In this system the coordinates of the Sun are 
(x,y,z)=(0. 0,0. 0,0.0) kpc and the centre of the Galaxy is 
located at (x,y,z) = (+8. 0,0. 0,0.0) kpc, having assumed that 
the distance to the Galactic center is 8.0 kpc. To avoid any 
possible effect of local inhomogeneities in the vicinity of the 
Galactic Plane we also removed the stars having \z\ < 0.5 
kpc. Given the adopted selections this is equivalent to ex- 
cluding stars within ~ 1 kpc from the Sun. 

From these samples we computed the density of RC 
stars on an x,y grid with boxes of size 1 kpex 1 kpc, spaced by 
0.5 kpc both in x and in y. The overlapping between adjacent 
boxes ensures some degree of smoothing, providing more 
clearly readable maps. In this way we obtained x,y density 
maps of the Southern and Northern Hemisphere that can be 
subtracted, "pixel-to-pixel", to obtain a map of the residual 
surface density in the x,y plane, that is residual maps of the 
Galactic Plane as seen from the Galactic pole. The density 
at each point of the subtracted map is computed as 



/ South North\ 

P(x,y) — \ n (x,y) - n (x,y) ) 



Note that nf^-f and n^°^ h ) are a sort of "column" density, 
since at the position (xo;yo) the derived density is equal to 
the number of RC stars having xq — 0.5 ^ x ^ xo + 0.5 
and yo — 0.5 ^ y ^ yo + 0.5, at any z within the limits 
of our sample. The imposed selection in latitude (5.0° < 
\b\ ^ 15.0°) imply that the range in z sampled by different 
pixels of the grid varies with their heliocentric distance. In 
particular, at Dq = 5.0 kpc the sampled range of z is 0.5 kpc 

~ \z\ ~ 1.5 kpc, at D Q = 10.0 kpc: 0.9 kpc ~ \z\ ~ 2.6 kpc, 
at Dq = 13.0 kpc: 1.1 kpc ~ \z\ ~ 3.3 kpc. It is important 
to realize that this has no consequence for the subtracted 
densities P( x , y ), since each pair of pixels to be subtracted 
have exactly the same volume. 

Since a large absolute density residual (e.g. in number) 
does not necessarily imply a statistically significant signal, 
it is mandatory to obtain the subtracted density maps also 
in terms of Signal-to-Noise ratio, that is in units of a: 



/ South North\ / 

P(x,y) — ( n (x, y) ~ n (x,y) )/ a (x,y) 



where 



T (x,y) 



(x,y) 



(x,y) 



Subtracted maps in terms of number density provide 
a physical idea of the extent of the detected overdensities, 



while subtracted maps in a units allow us to discriminate 
structures that are statistically reliable and significant from 
those that can arise from mere statistical fluctuations. 



3 OVERDENSITIES IN THE OUTER DISC 

Fig. 3 shows the maps of residuals obtained by subtracting 
the Northern hemisphere density map from the Southern 
hemisphere density map, in terms of number density (panel 
a) and in units of a (panel b). We recall here that the regions 
0° < I sC 40° and 320° «S I < 360° are excluded from the 
present analysis. In Appendix A it is shown that while the 
effect of contamination by dwarfs should not significantly 
affect subtracted density maps for K ^ 14.0, it may be- 
come quite relevant at fainter magnitudes. In this section 
we present subtracted maps that include also stars having 
14.0 < K < 15.0 to avoid abrupt distance thresholds in the 
proximity of the regions of interest. We caution the reader 
that the reliability of the maps may be drastically reduced 
beyond the heliocentric circle corresponding to the K — 14.0 
limit that is always clearly shown in Figs. 3, 4 and 5 below, 
as a reference. According to the results of Appendix A, the 
contamination by dwarfs can seriously affect the maps also 
around I = 180°. On the other hand, the main structures we 
deal with in the following analysis lie safely in the range of 
latitudes and distances that should be correctly probed by 
our subtracted-density technique (see Appendix A). 

With the only exception of a few limited low- 
significance structures near the edges of the map (in par- 
ticular near I = 40°, probably due to edge effects), the 
strongest South-North overdensities are confined in the third 
and fourth Galactic quadrants (180° ^ I ^ 360). This is 
not unexpected since it is well known that the Galactic 
Disc is warped toward the South Galactic hemisphere for 
180° ^ I ^ 360 an d toward the North Galactic hemisphere 
for 0°^ig 180 Jpiorgovski fc Sosinlll989l: lYusifovll2004l : 
lLopez-Corredoira et alJl2002l : IVig. Ghosh fc Oihal2005l and 
references therein) . The third and fourth Galactic quadrants 
appear dominated by a wide, galactic-scale overdensity that 
extends over the whole half of the disc with I > 180°. On top 
of this large-scale asymmetry there is a strong and spatially 
confined overdensity, with an elongated, approximately ellip- 
tical shape, with a denser roundish core located at (x,y)=(- 
3.7 kpc, -6.4 kpc), toward I = 240°, at R G c ~ 13.0 kpc and 
-Dq ~ 7.5 kpc. This is the Canis Major structure, as de- 
scribed in Pap-I, Pap-II and Pap-III. Note that a secondary 
core that appears towards I ~ 280° in Fig 3a is revealed as 
a non-statistically-significant structure by the inspection of 
the map in a units (Fig. 3b). 

The linear "finger of God" structures between I ~ 200° 
and I ~ 220°, and at I ~ 250° correspond to asymmetric 
features in the distribution of the extinction (see Fig. 2) 
and are likely artefacts associated with the (unknown) three 
dimensional distribution of the interstellar dust. These fea- 
tures are present in all the South-North maps presented in 
this paper and will not be discussed further. Note also that 
the strongest and most significant parts of the detected over- 
densities (including the body of CMa) lie within the circle 
enclosing the RC stars with K ^ 14.0 (inner heliocentric 
long-dashed circle, at Dq = 12.6 kpc), that is in the region 
where the contamination by dwarf stars is less severe (i.e., 
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Figure 3. Subtracted density maps (South-North) in terms of 
number density (panel a) and in units of a (Signal-to-Noise ratio; 
panel b). The direction of Galactic longitudes from I = 40° to 
/ = 320°, with steps of 20°, are indicated as labelled dotted lines. 
The encircled dot marks the position of the Sun, the X symbol 
indicates the Galactic Centre. Long-dashed circles are centred on 
the Sun: the outer circle indicates the distance limit of our anal- 
ysis, corresponding to RC stars having K = 15.0, the inner circle 
encloses the region sampled by RC stars having K $C 14.0, i.e. the 
range in which the contamination by dwarfs should be moderate. 
The continuous circle is centred on the Galactic Centre and has 
a radius Rqc = 13.1 kpc, corresponding to the heliocentric dis- 
tance of Dq =7.2 kpc for the centre of the Canis Major structure 
at I = 240° (Pap III). The region within ±40° in longitude from 
the Galactic Centre is excluded from the present analysis. 



Figure 4. Subtracted density maps (North-South) in terms of 
number density (panel a) and in units of a (Signal-to-Noise ratio; 
panel b). The symbols are the same as in Fig. 3. 



the c ontribution of dwarfs is less th an 50 per cent, accord- 
ing to lLopez-Corredoira et al . 2002, and the contribution of 
dwarfs to subtracted densities in lower than 30 per cent, see 
Appendix A) . 

Fig. 4 shows the maps of residuals obtained by subtract- 
ing the Southern hemisphere density map from the Northern 
hemisphere density map, in terms of number density (panel 
a) and in units of a (panel b). Also in this case the de- 
tected overdensities are confined to the Galactic quadrants 
that are expected to host the northern lobe of the Warp. 
The observed asymmetry has a large scale - as in the South- 
North maps - but is weaker and less significant than the 
large-scale Southern asymmetry. This may be due to dif- 
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ferences in the structu re of the Northern and the Sout hern 
Warp (see, for instance lYusifovl2004l : |p"orcel et al and 
to the effect of the disp lacement of the Sun w ith respect to 
the Galactic Plane (Ha mmerslev e t al. 1995)). However, in 
the present context is interesting to note that (a) there is no 
overdensity comparable to CMa in the Northern Hemisphere 
and (b) the most significant asymmetries lie in the range 
100° < I < 140°, that is at ~ 180° from the densest part of 
the large-scale Southern overdensity (280° ^ I ^ 320°). 

The inspection of Fig. 3 and Fig. 4 clearly demonstrates 
that CMa is the strongest (35 a) overdensity of the whole 
Galactic Disc. Moreover it is strongly spatially confined (the 
two innermost density contours have a typical scale of ~ 4 
kpc x 2 kpc) and there is no corresponding structure either 
on the opposite side of the disc (I ~ 60°) or at the reflected 
position about the y or x axis. 

Note that the detected excess of RC stars is quantita- 
tively very similar to what was found for M giants in Pap-I 
(in terms of a and of fractional excess). CMa appears as a 
strongly concentrated, cored overdensity on top of a galactic- 
scale asymmetry that has all the characteristics of the Galac- 
tic Warp. 

3.1 Removing the large scale asymmetry 

If the global smooth overdensity is related to the Galactic 
Disc (i.e. it reflects the Warp asymmetry) it should be pos- 
sible to subtract it from the maps of Fig. 3. Since the overall 
number of southern RC stars in our sample is 1.2 times that 
found in the Northern hemisphere, as the simplest assump- 
tion we rescale the density of the Northern map by this fac- 
tor and repeat the subtraction. The adopted normalization 
should account also for the effect of the off-plane position 
of the Sun. The resulting maps of residuals are shown in 
Fig. 5, both in terms of number density (panel a) and in cr 
units (panel b). Adopting a safe 10cr cut, the only surviving 
structure is CMa plus tiny traces of the already discussed 
finger-of-God features at / ~ 220°. 

Note that large negative residuals 10a) appear only 
near the edges of the maps (Fig. 5b). This means that the 
rescaled northern density map is a reasonable model of the 
large-scale overdensity observed in the southern hemisphere. 
We tried different simple models fitting the observed over- 
all decline of the large-scale overdensity with Rgc- in all 
cases the only overdensity surviving the subtraction of the 
model from the southern density map is CMa, in all cases 
displaying the characteristic shape shown in Fig. 5a, b. This 
exercise strongly supports our interpretation of the south- 
ern overdensities in the considered region of the Galaxy in 
terms of a large-scale Galactic asymmetry (warp) plus an 
additional unexpected elongated compact structure with a 
dense core, i.e. the Canis Major overdensity. 

3.2 Comparison with a Warp model 

In response to the suggestion that the CMa overdensity 
may be associated with the Galactic Warp, we have com- 
pared the observed characteristics of CMa with two dif- 
ferent models of the Galac tic Warp, namely t hat enclosed 
in the Galactic model bvjRobin^^jJ i2003T) . in Pap-II, 
and the warp model bv lYusifovri2004T^ in Pap-III. In both 
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Figure 5. Subtracted density maps (South-North) in terms of 
number density (panel a) and in units of a (panel b). In this case 
the densities in the Northern Hemisphere maps have been rescaled 
by a factor xl.2, to attempt to remove the large scale asymmetry 
that we identify as due to the Galactic Warp. In panel (b) we 
reported as small filled squares the pixels of the map where nega- 
tive residuals larger than 10c are detected. All the other symbols 
are the same as in Fig. 3. 

cases we found that the models were clearly unable to re- 
produce the very compact structure of CMa. In particular 
the distance profile of the overdensities predicted by mod- 
els are much more extended than what is observed for Ca- 
nis Major. Here we extend the comparison to the model by 
lLopez-Corredoira et alJ ll2002l hereafter L02) that has been 
constructed to fit the overall distribution of RC stars in the 
Galactic Disc. 

In order to compare density distributions obtained with 
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Figure 6. Subtracted density maps (South-North) from the observed sample (upper panel) and from the lLooez-Corredoira et all (2002) 
model (lower panel). The Model map has been normalized to have the same maximum density as the observed one. The inner long-dashed 
circle have a radius r = Dq = 7.2kpc, i.e. the heliocentric distance of CMa (Pap-III). The outer long-dashed circles encloses the region 
populated by RC stars having K $J 14.0. Both samples have also been limited to Rgc < 17.0 kpc. The continuous-line galactocentric 
circle plotted in the lower panel of the figure, having radius = 13.0 kpc, shows that the relevant comparison is performed inside the 
extrapolation region, hence the adopted form of the extrapolation cannot affect our results. Black dots in the lower panel indicates 
negative residuals larger than 50%, corresponding to the northern lobe of the Warp. 



the same technique, we constructed a Montecarlo realisation 
of the L02 model and we applied the same selections as to 
the observed sample (i.e., 5° |6| ^ 15° and zj ^ 0.5 kpc). 
We extended the L02 model outside its maximum Galacto- 
centric distance (Rgc ~ 13.0 kpc) with the following as- 
sumptions for the average elevation of the disc above the 
plane z w (see Eq. 20, in L02): 

z w = C w R[kpc] 5 - 25 sin($ - $ w ) + 15.0 



for 13.0 kpc < Rgc < 16.0 kpc, 
where Cw = 0.0012 and <3>vv = —5.0 deg, and 

z w = 0.0 

for Rgc > 16.0 kpc. This extrapolation produces a 
regime of nearly constant elevation for 13.0 kpc < Rgc ^ 
16.0 kpc and a drop to the plane at larger galactocentric dis- 
tances, that appear as the most appropriate for the Southern 
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warp, according IPorcel et alJ il997f> and lYusifovl (l2004f) . for 
example. For homogeneity, we applied the Rgc = 17.0 kpc 
cut also to the observed sample. Note however that the re- 
gion relevant for the comparison is largely included within 
the Rgc = 13.0 kpc limits, hence the adopted extension 
of the L02 model must have little effect on our results (see 
Fig. 6, below). 

Finally, we limit the comparison to stars within Dq = 
12.6 kpc, corresponding to the K ^ 14.0 for RC stars. After 
applying all of the above selection criteria, the Montecarlo 
realisation of the L02 model contains 10 6 particles. Consid- 
ering these particles as the analogues of RC stars, we have 
obtained subtracted density maps from the model in exactly 
the same way as we did for real stars (see Sect. 2.2). We 
normalized the density scale of the Model map by requiring 
that the maximum density is the same in the syntethic and 
observed map. We have tried also several other normaliza- 
tions and we found that the essence of the results remains 
unchanged. 

In Fig. 6 the observed South-North subtracted map (up- 
per panel) is compared to the same map obtained from the 
L02 model. Both maps are represented in the same scale of 
density excess. The differences between the two maps are 
clearly evident. In the model map the Warp appears as a 
large wedge-shaped overdensity, somewhat similar - in shape 
- to the densest part of the large scale overdensity observed 
in Fig. 3. The density of this feature reaches a broad max- 
imum in the range 270° ^ I ^ 280°, over a large range of 
Heliocentric distances, around (x,y)[ fcpc ] ~ (1.5;-10.3) - while 
the maximum of the observed density occurs at I ~ 240° and 
(x,y)[fc pc ] ~ (-3.5;-6.2), more than 6 kpc apart. It is quite 
clear that the model is totally unable to reproduce the CMa 
structure presented in the upper panel of Fig. 6 as the only 
significant overdensity, extending from I ~ 200° to I ~ 275° , 
with its compact elliptical shape and short minor axis, and 
with its core at / = 240°. The differences between the two 
structures (CMa and the L02 Warp model) clearly suggests 
that they are of a different nature and cannot be identified 
with one another, confirming the results presented in Pap-II 
and Pap-III (see also Fig. 1, here). 

Despite the large differences shown in Fig. 6, it may 
still be argued that the L02 Warp model still predicts some 
South-North overdensity at the position of CMa. Are the ob- 
served and predicted distributions locally similar, at this po- 
sition? To test this possibility we recur to the optical CMD 
of the XMM field (F-XMM, at I = 244.2° and b = -8.2°) we 
presented in Pap-II, that is reproduced in the middle panel 
of Fig. 7. To have an idea of how the CMD of an ensemble of 
old stars distributed according to the L02 model would ap- 
pear we proceeded as follows: (1) we take as template popu- 
lation the photometry of ~ 13000 stars of the globu lar clus- 
ter 47 Tucanae presented bv lBellazzini et all ll2004rj) . (2) we 
extracted, at random, 20000 stars from this sample and we 
corrected their magnitude for extinction and for distance 
(assu ming E(B-V)=0.04 and (m-M) =13.29 for the custer, 
from iFerraro et alJll99^ . (3) we assigned to each star a 
distance randomly extracted from the distance distribution 
predicted by the L02 model toward (l,b) = (244.2°, -8.2°), 
and (4) we corrected for the corresponding distance moduli 
and for the reddening of the F-XMM field (E(B-V)=0.14, 
see Pap-II). 

The synthetic CMD thus obtained is shown in the up- 



per panel of Fig. 7. Its morphology is completely differ- 
ent from that of the observed diagram. Stars around the 
Turn Off (TO) point of the cluster form a vertical sequence 
from V ~ 17.0 to V ~ 20.0 that abruptly bends to the 
red for V > 20.0 in the way typical of CMDs of Galactic 
fields that sample populations at widely different distances 
along the line of sight. This result is similar to what we ob- 
tained in Pap-II while comparing the F-XMM CMD with 
the predictions of t he warped and flared Galactic model by 
iRobin et ai] <2003l) . Note that ~30% of the stars plotted 
in the upper panel of Fig. 7 have associated distances that 
place them in the regime in which we adopt an extrapola- 
tion of the L02 model. However, repeating the experiment 
while forcing the selection of stars with Rgc ^ 13.0 leads 
to a CMD having a morphology essentially indistinguishable 
from that shown here. Hence, also in this case the adopted 
extrapolation of the L02 model has a negligible effect on 
our results. In contrast, the F-XMM CMD shows a clear 
and narrow Main Sequence that bends continuously from 
the TO (around V ~ 18.5) to the limit of the survey, clearly 
resembling a bound stellar system. This similarity is made 
clearer by the comparison with the CMD shown in the lower 
panel of Fig. 7. This synthetic diagram has been obtained 
repeating the steps 1, 2 and 4 described above but assign- 
ing to stars distances extracted from a Gaussian distribu- 
tion having a mean Dq = 7.2 kpc (i.e., the distance of 
C Ma as derived in Pap-III) and a = 0.8 kpc, as derived 
bv lMartinez-Dekado et alJ ll2005al) for CMa. Except for the 
foreground field population that is, obviously, present in the 
F-XMM and absent in the higher latitude 47 Tuc field, the 
two MS features are strikingly similar. This test confirms 
that the observed CMa overdensity in not predicted by cur- 
rent models of the Galactic Warp, neither on global or lo- 
cal scales. While all the considered models predict that the 
Galactic Warp should provide some overdensity at the po- 
sition of CMa, in good agreement with the results shown in 
Fig. 3, none of them is able to reproduce anything remotely 
similar to the extremely compact, cored structure we observe 
in Canis Major. 

It should be recalled that the Galaxy models including 
a warped and flared outer Disc, as those considered here and 
in Pap-II and Pap-III, are not meant to provide a perfect 
representation of the Galaxy, in particular for what con- 
cerns the outer parts of the Galactic Disc where the avail- 
able observations are scanty and uncertain over large ranges 
of longitude. Hence, the fact that these models are unable 
to reproduce the observed structure of CMa is not an ulti- 
mate proof that Canis Major cannot be a substructure of 
the Galactic Disc. On the other hand, this is a proof that 

(a) Canis Major is a truly unexpected structure, since no 
similarly spatially confined and dense structure is included 
in existing warped and flared models of Galactic Disc and, 

(b) the CMa structure is fairly different from the general 
expectations of an overdensity due to the Galactic Warp, 
since it is a compact elongated and cored structure instead 
of a large wedge-shaped lobe, as shown in Fig. 6, and, above 
all, it appears at a significantly different position. 

3.2.1 The Argo structure 

While the present paper was near completion a preprint 
was posted iRocha- Pinto et all2005l hereafter RP05) study- 
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Figure 7. Synthetic (upper and lower panels) and observed (middle panel) Colour Magnitude Diagrams towar d (l\b) = (244.2° ; —8.2°) . 
Upper panel: CMD of ~ 20000 stars randomly extracted from the photometry of the cluster 47 Tucanae by iBellazzini et~ . 12004b), 
as it would appear if the considered stars were distributed according to the distance distribution predicted by the Lopez-Corredoira at 
al. Warp model in the above direction. Middle panel: observed CMD toward (l;b) = (244.2°; -8.2 °) (XMM Field, f rom Pa p-II). Lower 
panel: CMD of ~ 20000 stars randomly extracted from the photometry of the cluster 47 Tucanae by Bellazzin i et all l2004bT) . as it would 
appear if the considered stars were distributed according to a Gaussian distrib ution having mean Dm = 7.2 kp c and standard deviation 
a = 0.8 kpc, i.e. the line of sight distance distribution of CMa as estimated by Martincz-Dclg ado et alj J2005al) . 



ing the distribution of 2MASS M giants at low galactic 
latitudes (\b\ ^ 20°). Adopting an extinction grid with a 
finer resolution compared to previous an alyses of M-giants 
ijMaiewski et alJl2003L ICrane et alJl2003t Pap-I), and using 
a starcounts model to subtract the smooth Galactic com- 
ponents, they are able to follow overdensities down to lower 
absolute latitudes (\b\) than, for instance, in Pap-I. They find 
that, in the 220° ^ / ^ 320° region, all the existing over- 
densities are enclosed in the range of Heliocentric distances 



6.0 kpc < Dq < 20 kpc. They confirm the presence of the 
CMa overdensity in this range but they claim that it is part 
of a larger and denser structure, centred around I ~ 290°, 
that they interpret as a large dwarf galaxy which they name 
Argo. They suggest that Argo and CMa can be erroneously 
seen as separated entities only because of a plume of extinc- 
tion located at / ~ 265° extending to b ~ —12° (see Fig. 2 
and Fig. 8, below) that introduces an artificial discontinuity 
in the surface density. 
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In Fig. 8 we present the South-North subtracted surface 
density map of RC stars in the range 200° ^ I ^ 320°, 
obtained by counting stars (in the distance range 5.0 kpc 
~ Dq ~ 13.0 kpc) within 4° x 4° boxes spaced by 2° 
both in longitude and latitude. In this case also stars with 
\b\ < 5.0° have been included in the sample, for a better 
comparison with RP05. The main body of the claimed Argo 
structure is clearly traced, in surface density, also by RC 
stars. In Fig. 8 it appears as a strong, elongated overdensity 
with a broad maximum around / ~ 290° — 310°, with CMa 
placed at its low longitude edge (compare with Fig. 1 and 
Fig. 3 of RP05). 

The maps displayed in Fig. 3a,b (as well as Fig. 5, 6, 
and 7, above) show that there is no spatially confined 3-D 
structure corresponding to this strong surface overdensity. 
Integrating over the range 6.0 kpc < Dq < 20 kpc along the 
I ~ 280° — 300° direction implies summing up the contribu- 
tion of all the stars encountered along the line of sight that 
crosses the maximum of the overdensity produced by the 
Galactc warp. In any case, while in terms of surface density 
- integrated over 14 kpc along the line of sight - CMa is a 
much weaker structure compared to Argo (see Fig. 8) , it cor- 
responds to a much more spatially confined structure: CMa 
spans the range 5 kpc < D aun < 12 kpc, with a line-of-sight 
Half Width at H alf Maxi mum o f just ~ 2.0 kpc (see Pap-II, 
iMartinez-Deleado et"afl J2005al) and Sect. 4.1, below). The 
major axis of CMa is approximately oriented along the path 
of a nearly-circular orbit around the centre of the Galaxy 
(see Fig. 5 and Fig. 11, below), consistent with the N-body 
simulation presented in Pap- IV (see Sect. 4.2, below). On 
the other hand, the Argo structure corresponds to a part of 
the large-scale southern overdensity that spans more than 15 
kpc in heliocentric distance, toward / ~ 290° and it does not 
appear as a coherent, cored structure as CMa (see, Fig. 3). 

The above arguments are clearly displayed in Fig. 9, 
where we consider the surface density of RC stars confined 
to the very narrow range of Galactocentric distances includ- 
ing most of the main body of CMa (12.0 kpc ^ Rgc ^ 
14.0 kpc) . The subtracted surface density map in the upper 
panel of Fig. 9 has been obtained in the same way as Fig. 8, 
but is expressed in terms of fractional density excess to al- 
low a comparison with the L02 model. It is quite evident 
that at this spatial location the only strong overdensity is 
CMa, while Argo has completely disappeared. On the other 
hand the fractional excess due to the Warp (integrated over 
5° ^ |6| ^ 15°, Fig. 9, lower panel), as predicted by the L02 
model, reaches its maximum at the longitude corresponding 
to the maximum density of Argo (I ~ 290° — 300°). 

There is also a quantitative self-consistency problem to 
consider. RP05 derive (from M giants) an estimate of the 
total luminosity of L = 3.7— 15 x 10 6 Lq for Argo. Limiting 
ourselves to the stars contributing to the map of Fig. 8 and 
within 265° ^ I ^ 320° we count an excess of more than 
1.8 x 10 6 RC stars. With the assumptions stated in Sect. 4.3, 
below, we obtain a total luminosity of L ~ 1.2 x 10 9 , orders 
of magnitude larger than the R05 estimate. Note that our 
estimate should be considered as a lower limit since (a) the 
contribution of CMa is not included, while RP05 consider 
it as a part of Argo, and (b) we limited the calculation to 
the approximate distance range 5 kpc ^ D 3un 13 kpc, 
while R05 integrate over 6.0 kpc < Dq < 20 kpc. Hence 
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Figure 8. Subtracted map (South-North) in terms of surface den- 
sity in the sky for RC stars in the range 11.5 $C K $C 14.0. The 
surface density has been computed on a grid of 4° x 4° pixels 
spaced by 2° both in latitude and longitude. In this case stars 
with |b| < 5.0° have also been included in the sample. The posi- 
tion of CMa is indicated by a labelled arrow. Small squares mark 
positions where the reddening is 0.5 ^ E(B — V) < 1.0, large 
squares correspond to reddening in excess of 1.0 mag. 



the beautiful consistency among the luminosity estimates 
obtained from M-giants, RC stars and MS stars achieved 
for CMa (see Sect. 4.3, below) is completely lost in the case 
of Argo. This strongly suggests that Argo is not a coherent 
stellar system, but is instead just a surface density structure 
that originates from the mix of different (and unrelated) 
populations that happen to lie along the same line of sight. 

RP05 discard the warp hypothesis for the nature of 
Argo claiming that the shape of the structure is not "warp- 
like". However, in their view, all the southern overdensi- 
ties in the 220° ^ I ^ 320° region are due to Argo and 
CMa. The implicit conclusion is that the Galactic Warp 
(and flare) does not produce any observable South-North 
asymmetry in this region, at o dds with the conclu s ions o f 
lLopez-Corredoira et alJ j2002D : IVig. Ghosh fc Oihal J2005D . 
The smooth nature of the global South-North asymmetry 
shown to be present in this region, its Galactic-scale dimen- 
sions, and the fact that a rescaled version of Northern den- 
sity map is a good model for it (Fig. 5), strongly indicate 
that the surface density structure centred around I = 290° 
is not a dwarf galaxy remnant (Argo) but the effect of the 
Warp (and flare) asymmetries of the Galactic disc. 

Turning back our attention to CMa, there is an inter- 
esting point that can be appreciated from the comparison of 
Figs. 9 and 10. While the CMa structure is present in both 
maps, its appearance is quite different. Most of the difference 
is just apparent and is due to the density scale adopted in 
Fig. 8, that is best suited to provide a good view of "Argo" 
and compress the whole structure of CMa to just three lev- 
els of density. However part of the difference is real, and 
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Figure 9. Upper panel: subtracted map (South-North) in terms 
of fractional excess of surface density in the sky for RC stars in the 
range of Galactocentric distances 12.0 kpc ^ Rgc Ss 14.0 kpc, i.e. 
the range enclosing the densest part of the CMa structure. Lower 
panel: fractional density excess as a function of Galactic longitude 
in the same range of Galactocentric distances, as predicted by the 
Lopez-Corredoira et al. Warp model. Data for |6| < 5.0° are not 
included here. 

is due to the much larger degree of foreground/background 
contamination affecting the map of Fig. 8, which includes 
stars at any Galactocentric distance. In particular, the cen- 
ter of the structure (i.e., the region of maximum density) is 
shifted from I ~ 242° in the distance limited map (Fig. 9, 
that should be closer to the actual physical structure of the 
system) to I ~ 244° in the all-distances map of Fig. 8, due 
to contamination by likely unrelated stars at various dis- 
tances along the line of sight. This kind of effect should be 
always taken into account in studying the structure of CMa, 
since effects due to the integration along the line of sight can 
significantly change the overall view of the system. 

3.3 The nature of the Canis Major overdensity 

In our view there are only two possible explanations for the 
observational evidence presented above: 

(i) Canis Major is a strong and unexpected substructure 
of the outer Galactic Disc. This would imply that the Disc 
hosts very dense and compact substructures of old stars (the 
age should be larger than 1-2 Gyr to have RC stars and 
should be large r than 4-5 Gyr according the op tical CMDs, 
see Pap-II and iMartinez-Delgado et alJ i2005al) ). Such sub- 
structure should be very different, in nature, from usual spi- 
ral arms, for instance, that host populations dominated by 
very young stars. A local distorsion of the outer Disc, cou- 
pled with a large scale warp, may also possibly reproduce 
the observed distribution of overdensities in the Southern 
Galactic Hemisphere. 

(ii) Canis Major is the remnant of a disrupting dwarf 



galaxy, as suggested in Pap-I, II, III, and IV. This interpre- 
tation naturally fits all the observational evidence collected 
up to now, with particular emphasis on the size and den- 
sity of the structure and to the kinematics o f its sta rs (see 
Pap-III, Pap-IV and lMartinez-Dekado et alJ J2005bl) ). 

For the reasons described above and in the previous 
papers of this series, our preference goes to the second hy- 
pothesis, that appears to provide a much easier explanation 
of the whole observational scenario. In the following we will 
take this hypothesis as true and, on this basis, we will derive 
some interesting physical characteristics of CMa as a stellar 
system. 



4 PHYSICAL PARAMETERS OF CANIS 
MAJOR 

4.1 Distance 

In Pap-III we obtained a distance profile of CMa using M- 
giants as distance indicators. Here we can repeat that anal- 
ysis using the more reliable distance scale provided by RC 
stars. In the upper panel of Fig. 10 we show the distance 
distribution of the RC stars in our sample in the Northern 
(blue) and Southern (red) hemispheres, in a wide region to- 
ward the centre of CMa. The southern profile shows an over- 
all excess at any distance, with respect to the northern one, 
plus an additional bump peaked around Dq ~ 7.5 kpc. We 
rescaled the northern distribution by the same factor as in 
Fig. 5 (e.g., xl.2) and subtract it from the southern profile. 
The result is shown in the lower panel of Fig. 10: the South- 
North excess appear to be confined within Dq — 5 kpc and 
-Dq ~ 12 kpc. The peak of the profile is at Dq = 7.2 ±1.0 
kpc and the Half Width at Half Maximum is HWHM ~ 2.0 
kpc, in excellent agreement with the re sults of Pap-I, Pap- 
III and Martinez-Delgado et alJ (1200 5a). There are a few in- 
teresting considerations that can be deduced by inspecting 
Fig. 10: 

• The distance profile of CMa as obtained from its RC 
stars is strikingly similar to that derived from M-giants in 
Pap-III. Even leaving aside the perfect coincidence of the 
peaks of the two distributions, that may be ascribed to a 
mere coincidence within the uncertainties associated with 
the two distance scales, it remains a fact that the two inde- 
pendent tracers give essentially the same line-of-sight profile 
of the CMa system. 

• In our various analyses we have obtained distance esti- 
mates for CMa from four different and independent distance 
indicators. In order of accuracy: from the tip of the Red Gi- 
ant Branch in the I band (Dq ~ 7.2 kpc, Pap-III), from 
the Red Clump in K band (-Dq — 7.2 kpc, present anal- 
ysis), from the photometric parallax of M giants in the K 
band (-Dq ~ 7.1 kpc, Pap-I), and from Main Sequence fit- 
ting in the B and V bands (D ~ 7.6 - 8.7, Pap-II). All of 
these independent methods consistently provide a distance 
of the centre of CMa within 7.0 kpc and 9.0 kpc. Hence, the 
average distance of CMa should be now regarded as a very 
reliable piece of information about this system. Note that 
if Dq ~ 7.2 kpc is assumed, the age estimate obtained in 
Pap-II would shift toward older ages (i.e. 6 Gyr< age < 10 
Gyr, instead of 4 Gyr< age < 10 Gyr). 
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Figure 10. Distance profile of the Canis Major overdensity. 
Upper panel: distance distribution of RC stars in the range 
230° < I < 250° and 5° < |6| < 15°, in the Southern (Red) 
and Northern (blue) Galactic Hemisphere. Lower panel: residual 
of the subtraction of the rescaled Northern profile (by a factor 
xl.2) from the Southern profile. The long-dashed lines approx- 
imately encloses the range of variation of the flat part of the 
residual in response to variations of the scaling factor by ± 10 
per cent. 



4-1.1 The spatial orientation of CMa 

The subtracted density maps of Figs. 3, 5 and 6 make evi- 
dent that the physical orientation of the elongated body of 
CMa implies a heliocentric distance-longitude gradient: the 
heliocentric distance of CMa grows with Galactic longitude. 
Using the same technique as above we find that the dis- 
tance of the main body of CMa varies from Dq ~ 6.3 kpc 
at I ~ 225°, to D Q ~ 7.2 kpc at I ~ 240°, to D Q ~ 9.3 kpc 
at I ~ 265°. 

The effect is clearly depicted in Fig. 11, where the con- 
tours of the main body of CMa (in red, taken from Fig. 5b) 
are superposed on the results of the N-body simulation of 
the disruption of the CMa galaxy presented in Pap-IV. The 
simulation is produced by a dwarf galaxy model of mass 
5 x 10 8 Mq that is accreted onto the Milky Way during ~ 3 
Gyr (see Pap-IV for details) and it was built to reproduce 
the current mean position and 3-D velocity of CMa as well 
as the distance-radial velocity gradient discovered in Pap- 
IV. Hence the distance-longitude behaviour of the remnant 
is a prediction of the model, not an a priori requirement. 
Therefore the excellent qualitative and quantitative agree- 
ment with the observed and predicted spatial orientation of 
the main body of CMa may be considered as a success of 
the dwarf-galaxy-accretion scenario for the nature of CMa. 
In other words, the structure and spatial orientation of CMa 
as obtained from RC stars is fully consistent with the model 
of a disrupting dwarf galaxy that fits the observed position 
and kinematics of CMa. 



Figure 11. The excess-density contours of CMa taken from 
Fig. 5b (down to 16a, in red) are compared with the N-body 
model of the disruption of the CMa dwarf galaxy presented in 
Pap-IV (small dots). 



4.2 Luminosity and surface brightness 

The Evolutionary Flux relation iRenzini fc Buzzonil Il986l : 
iRenzini fc Fusi Peccill988t lRenzinill998l) states that the to- 
tal number of stars in a given evolutionary phase j (Nj) is 
proportional to the total luminosity of a stellar system (Lt). 
Considering a Simple Stellar Population (SSP, i.e. a popula- 
tion of stars with the same age and chemical composition): 



N 3 = B(t)L T tj 



(4) 



where tj is the duration of the evolutionary phase in years, 
and B(t) is the specific evolutionary flux, i.e. the number of 
stars entering or leaving any post Main Sequence evolution- 
ary phase per year per solar luminosity as a function of time 
(t, i.e. the age of the SSP). 

We can use Eqn. 4 to obtain a distance-independent 
estimate of the total luminosity of CMa from the RC star 
counts. With some algebra we obtain: 



M v = -2.5 



log[i^]-log[^]-logfe]-log[B(t)] 
Lv 



+M vQ (5) 



where My is the absolute integrated magnitude in the V- 
band and M v q is the absolute V magnitude of the Sun. 
For the present case we adopt B( t) = 1.71 x 10~ n and 
-§21 = 1.5 from the SSP models by iMarastor] Jl99cl 120051) 
with [M/H] = -0.33, age= 6.0 Gyr and Salpeter Initial 
M ass Function 1 . M v q = 4.83 is also adopted according 
to iMarastonl 1^998:), while the duration of the core He- 
burning phase (tuB = 0.87 x 10 8 yrs) for a population with 
age~ 6.0 Gyr and \MI H\ = —0.35 has been taken from 
iPietrinferni et alJ (120041) 2 . 



1 see http://www.mpe. mpg.de~maraston/SSP 
see |http://www.te.astro.it/BASTI/| " 
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In the same way we can use Eqn. 5 bv lRenzinl lll998t) to 
obtain a distance-dependent estimate of the surface bright- 
ness near the centre of CMa (/j,v,o)'- 

Hv.o = -2.5 log(ry) + 2.5 log(FoV) + (to - M) + Norm(6) 

where rij is the number of stars in the j evolutionary phase 
within the considered Field of View (FoV, in arcsec 2 , for ex- 
ample) , and the term Norm encompasses all the theoretical 
assumptions: 

Norm = 2.5 log[^-] ± 2.5 logfo] + 2.5 log[B(t)] + M v& 
Ly 

We assume the distance modulus (m — M)o = 14.3, as de- 
rived in Pap-III and in the present contribution. 

To estimate the observables required by Eqn. 5 and 
Eqn. 6 (Nhb and uhb) we proceeded as follows. We se- 
lected from our sample the RC stars of both hemispheres 
within 200° < I < 280°, 15.0° < \b\ < 5.0° and 10.0 kpc 
^ Rgc ^ 16.0 kpc (the adopted selection is displayed in 
Fig. 12). Then we subtracted the total number of North- 
ern RC stars - rescaled by the usual 1.2 factor - from the 
total number of Southern RC stars in this region. We ob- 
tain an excess of ~ 1.1 x 10 5 RC stars that can be as- 
sociated to CMa, corresponding to a total luminosity of 
Lt — 7.4 x 10 7 Lq. The integrated absolute magnitude is 
My = —14.4 ± 0.8, where the uncertainty includes a ±10 
per cent variation in the xl.2 scaling factor adopted to take 
into account the contribution of the large-scale overdensity 
(Warp) to the detected excess. This estimate is in excellent 
ag reement with those obtai n ed in P ap-I from M giants and 
by iMartinez-Deleado et al . (2005a) from optical photome- 
try (My = —14.5 ±0.1). We estimate that the inclusion or 
exclusion from the sample of the stars within the |6| < 5° 
region may change the value of the V-band integrated mag- 
nitude by ~ ±0.5 mag. 

To estimate the surface brightness we computed the 
South-1.2x North excess of RC stars within a 2° x 2° box 
centred on (l,b) = (240°, -6°) and on (l,b) = (244°, -6°). 
We chose two different positions to provide estimates of the 
surface brightness for the centre of the structure as obtained 
from a distance-selected sample (as in Fig. 9) and from the 
straight surface density distribution, as in Fig. 8. We obtain 
a density of ~ 314 RC stars deg -2 at I = 240° , corresponding 
to fiy = 24.0 ± 0.6 mag/arcsec 2 , and ~ 370 RC stars deg -2 
at I — 244°, corresponding to /iy = 23.8 ± 0.6 mag/arcsec 2 . 
T hese values are in reason a ble agr eement with the results 
of IMartinez-Deleado et al ] J2005ah obtained from data at 
{I, b) = (240°, -8°) Ov.o = 23.3 ± 0.1 mag/arcsec 2 ). 

Note that B(t), at fixed [M/H] = -0.33, varies from 
1.66 x 10" 11 at t = 1.0 Gyr to 1.75 x 10 -11 at t = 13.0 Gyr, 
and, for fixed age = 6.0 Gyr, it ranges from 1.66 x 10 -11 
at [M/H] = -1.35 to 1.98 x 10~ n at [M/H] = ±0.35. The 
variation of Ihb over th e same range of age and metallicity 
is just of a few per cent JPietrinferni et al.ll2004lb Hence the 
particular choices of the above parameters do not seriously 
affect our final estimates. Other factors, not completely ac- 
counted for, as the uncertainty in the exact position of the 
center of CMa, the overall extent of the system, and the 
adopted SSP hypothesis, are expected to be the dominant 
contributors to the uncertainty of our estimates of the in- 
tegrated magnitude and of the central surface brightness of 
the system. 
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Figure 12. Subtracted density map taken from Fig. 5a (North- 
ern Hemisphere densities rescaled by a factor xl.2). The contin- 
uous circles with radius -Rgc = 10-0 kpc and Rgc = 16-0 kpc, 
and the continuous lines indicating the directions I = 200° and 
I = 280° enclose the region that we have considered for the com- 
putation of the total number and the central surface density of 
RC stars attributable to Canis Major. These numbers have been 
used to compute the integrated absolute magnitude and the sur- 
face brightness (toward two different directions) of the system, 
that are reported in the figure. 

CMa turns out to have a total luminosity and a surface 
brightness quite similar to those of the Sgr dSph, as already 
noted in Pap-I, but the above results suggest that it may be 
slightly more luminous and with a brighter fly, 0, in agree- 
ment with t he My vs. fJ,y,0 relation of d warf spheroidal 
galaxies (see iMartinez-Delgado et al.ll20 05a. and references 
therein). 

It is interesting to note the striking effect of the relative 
distance on the number densitie s for two stellar syste ms of 
similar central surface brightness. iMonaco et all (120031) mea- 
sured a density of RC stars of ~ 2100 deg -2 in the central 
degree of the Sgr dSph, while here we observe just 314 deg -2 
at / = 240° , yet the surface brightness of the two systems is 
similar. This may largely account for the great difficulty in 
the identification of the RC in optical CMDs from off-centred 
fie lds of 0.25 deg 2 or less as those presented in Pap-II and 
in iMartinez-Delgado et all (|2005a). This argument will be 
developed further in Sect. 4.4, below. 



4.3 How should a nearby galaxy appear? 

CMa is by far the nearest dwarf galaxy known, a factor 
~ 3.6 closer than the previous "record holder" , the Sgr dSph 
(Mon aco et alJl2004fl . This is the main reason why it ap- 
pears so unusual, at first glance: a structure covering tens 
of degrees on the sky and with as few as ~ 300 — 400 RC 
stars in its central square degree. This is very far from our 
idea of the typical dwarf satellite of the Milky Way, that 
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has been acquired by studying much more distant exam- 
ples. To show that there is nothing unusual in CMa except 
its distance we undertook the exercise of seeing how the Sgr 
dSph galaxy would appear if it were placed at the same dis- 
tance and position as CMa. In the following we will adopt 
the structural parameters of Sgr (core radius r c , limiting ra- 
di us rt, central surface br ightness and axis ra tio) as derived 
bvlMaiewski et all |2008), the distance from iMonaco et alJ 
(2004) and the sur f ace d ensity of Horizontal Branch stars 
from lMonaco et alJ J2003h . 

Fig. 13 shows the core and the limiting contours of Sgr 
once placed at the distance and position of CMa. The up- 
per panel of Fig. 13 has the same longitude scale as Fig. 9, 
to allow a direct comparison. The core of the "mock" Sgr 
shows a major axis of more than ~ 26° (r c — 13.3°) and 
the overall appearance is strikingly similar to CMa. Note 
that the limiting contour extends well beyond I = 200° on 
the low longitude side, and beyond I = 300° on the other 
(r t = 64.5°). Note also that the outskirts of the "mock" Sgr 
galaxy extends into the Northern Galactic Hemisphere, up 
to 6- +5/10°. 

From Eqn. 6 the following relation can be derived be- 
tween the surface brightness (at any chosen wavelength A) 
A*a,2 and the number of stars per given area in a given 
evolutionary phase (ni,n2) for two stellar systems located 
at distances d\ and d^: 

- Ma, 2 = -2.51og(rai/n 2 ) + 51og(di/d 2 ) (7) 

The second term on the right hand side of Eqn. 7 dis- 
plays the effect of distance on the number surface density of 
a system of a given surface brightness 3 . To rescale the num- 
ber surface density of Sgr from di — 26.3 kpc to d\ — 7.2 kpc 
we must assume /ix,i — Ma, 2 = 0.0. In the present case we find 
that the observed surface number density of Sgr stars should 
be reduced by a factor 0.075, if the galaxy is placed at the 
distance of CMa. Hence, in this case, we would observe ~ 150 
RC stars deg -2 , 10 RR Lyrae deg -2 and 12 Blue Horizon- 
tal Branch stars deg -2 in the centre of the Sgr galaxy. This 
effect fully accounts for the difficulty of detecting clearly- 
defined sequences of evolved stars in small fields (smaller 
than 1 deg 2 ) even within the main body of CMa, as well as 
for the low yield in searches of rare stellar species (as RR 
Lyrae stars in a system domin ated by an intermediate-age 
and metal rich po pulation, see lKinman. Saha fc PieHl2004l; 
iMateu et ai]l2005l) . 

In conclusion, the apparent size of CMa and the number 
density of its stellar populations are fully compatible with 
those of a dwarf galaxy similar to Sgr and located at 7.2 kpc 
from us. 



4.4 Old open clusters in CMa 

The possible association of some old open cluste rs with CMa 
was already suggested in Pap - 1 (see also Pap-II: ICrane et alJ 
120031 : iFrinchabov et alJ l2004). A thorough analysis of these 

3 It has to be noted that the effect of the physical size of a sys- 
tem on the spread of sequences on a CMD weakens with increas- 
ing distance. This is an additional factor that favours the easier 
recognition of distant stellar structures compared to more nearby 
systems. 
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Figure 13. The appearance of the Sgr dSph if it were placed 
at the distance and position of CMa. The inner ellipse encloses 
the core of the galaxy (densely-shaded region), the outer ellipse is 
the limitin g contour ( sparsely shaded region; Sgr parameters from 
iMaiewski et alj l2003i. The upper panel has the same longitude 
scale as Fig. 8, to allow a direct comparison with the map shown 
there. 

clusters is clearly beyond the scope of the present contribu- 
tion. However, the use of the Red Clump as a tracer may 
provide the unique opportunity of placing a few clusters onto 
the same distance scale as the field population of CMa. In 
Table 1 we list a number of open clusters selected to be 
older than ~ 0.8 Gyr and to lie in the surroundings of the 
main body of C Ma (according to the WEBDA database, 
lMermilliodlll995Tl . In particular the first group presented in 
Tab. 1 is composed of three clusters that directly project 
onto the main body of CMa (as seen in Fig. 8 and 10, for 
instance), and for which it has been possible to estimate the 
K magnitude of the Red Clump from 2MASS photometry. 
Hence the distance estimates for To 2, AM 2 and Haf 11 
have been derived in exactly the same way and under the 
same assumptions (luminosity of the standard candle and 
reddening) used for the analysis of the field population per- 
formed above. These estimates strongly indicate that To 2, 
AM 2 and Haf 11 are physically located within the main 
body of CMa (see Pap-II for further discussion and refer- 
ences). Note that, while AM 2 and Haf 11 lack any radial 
velocity estimate, To 2 has V r ~ +114 km/s, fully consistent 
with the systemic velocity of CMa (Pap III, Pap IV). 



5 SUMMARY AND CONCLUSIONS 

We have performed a differential analysis of the distribu- 
tion of the Southern Galactic hemisphere excess of colour- 
selected Red Clump stars extracted from 2MASS in the re- 
gion 40° ^ I ^ 320° and 5° |6| ^ 15°. Our main aim was 
to trace the structure of the main body of the newly discov- 
ered Canis Major stellar system (Pap-I) and of the stellar 
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Table 1. Old open clusters around the main body of CMa. 



Name 


1° 


b° 


E(B-V) 


k rc 


(m - 


-M) 


D sun [kpc] 


Tombaugh 2 


232.83 


-6.88 


0.29 


12.95 ±0.10 


14.45 


±0.22 


7.76 ±0.8 


Haffner 11 


242.39 


-3.54 


0.57 


12.20 ±0.20 


13.70 


±0.28 


5.49 ± 0.8 


Arp-Madore 2 


248.12 


-5.88 


0.51 


13.65 ±0.10 


15.15 


±0.22 


10.71 ± 1.1 


Berekeley 36 


227.38 


-0.59 










6.1 


NGC 2243 


239.48 


-18.01 










1.1 


Melotte 66 


259.56 


-14.24 










4.3 


Ruprecht 75 


276.79 


-4.48 










4.3 


CC06 


238.48 


-4.28 










4.5 


Van den Bergh-Hagen 66 


276.00 


-3.01 










7.0 


Saurer 2 


257.99 


-1.01 










4.8 



a 



a All the reported data, with the except ion of those relati ve to the RC distance of To 2, AM 2 and Haf 11, 
are drawn from the WEBDA data base (Mcrmilliod 1995), sec http://obswww.unige.ch/webda/ Data for 
CC06 are from ilvanov et al.ll2002h . First group: clusters within 230° ^ I ^ 260°, -12° < b < 3° and 4.0 
kpc < D sun ^11.0 kpc, with age ^ 1 Gyr. Second group: clusters within 220° < / ^ 280°, -20° < b < 0° 
and 4.0 kpc sC D aun ^11.0 kpc, with age > 0.8 Gyr. 



component of the southern part of the Galactic Warp (Pap- 
II, Pap-III). 

We have obtained detailed maps of the South-North 
and North-South overdensities in the Galactic Disc. These 
maps show that the third and fourth Galactic quadrants 
(180° ^ / ^ 360°) are dominated by a smooth large-scale 
South-North asymmetry that extends all over the considered 
half of the Disc and over a large range of heliocentric dis- 
tances. The above characteristics, the comparisons with the 
large-scale North-South asymmetries observed in the first 
and second quadrants (0° ^ / ^ 180°) and the comparison 
with the warped and flared Galactic model of L02 strongly 
suggest that such a wide South-North overdensity is due to 
an asymmetry of Galactic scale, in particular the southern 
lobe of the Warp of the Galactic Disc. 

The Canis Major structure, as recognized in Pap-I and 
Pap-2, appears - in the South-North density maps - as the 
strongest spatial overdensity of the whole Galactic Disc, 
either in terms of number density or of statistical signifi- 
cance. Simply subtracting a rescaled version of the Northern 
hemisphere density map (by a factor 1.2) from the South- 
ern hemisphere density map, the large scale asymmetry at- 
tributable to the Galactic Warp is completely removed and 
Canis Major remains as the only significant overdensity, with 
a peak at more than 25 a. The structure appears elongated 
along the tangential direction, extending from I ~ 200° to 
/ ~ 280°. It has a nearly elliptical shape and it is strongly 
spatially confined, with a HWHM along the line of sight at 
its center of ~ 2.0 kpc (in excellent agreement with the re - 
sults of Pap-I, Pap-III, and lMartinez-Deleado et al.ll2005alh 
All the maps show that Canis Major has a roundish core, 
located at Dq = 7.2 ± 1.0 kpc, toward I = 241.7°, at a 
Galactocentric distance of Rgc ~ 13.1 ± 1.0 kpc. 

The L02 model of the Galactic Warp is unable to pre- 
dict the existence of Canis Major as a structure of Galactic 
origin, both on the large and on the lo cal scale . The same 
has been shown t o be true also for the lYusifovl l)2004t) and 
the iRobin et al ] (120031) models, in Pap-III and in Pap-II, 
respectively (see also Fig. 1, above). Within the framework 



set by the observational evidence listed above, Canis Major 
appears as a compact and confined stellar system that is 
"superposed" on the Warp structure (see Fig. 3, for exam- 
pie). 

We also investigate the c laim made in the recent study 
bv iRocha- Pinto et alJ (Ho05), that Canis Major is an ex- 
ternal field of a larger Argo structure, which would be the 
true centre of the dwarf galaxy remnant accreted onto the 
Milky Way. The distribution of red clump stars along the line 
of sight does not support this conclusion, however. Unlike 
the situation in the Canis Major region, we find no signifi- 
cant spatially-confined structure in the line of sight through 
Argo. This line of sight, instead, crosses the region of the disc 
were the large-scale overdensity associated with the Galac- 
tic Warp reaches its maximum. We conclude that Argo is 
most likely not related to Canis Major, and is probably an 
asymmetry of Galactic origin. 

The above observational scenario does not provide the 
ultimate word on the nature of CMa. Given our poor knowl- 
edge of the outermost region of the Galactic Disc and the 
obvious limitations in current models of the Disc itself, the 
possibility that CMa is a substructure of the Galactic Disc 
cannot be excluded yet. However, in this case, we should 
admit the presence of an unexpected large and very dense 
substructure and/or local distorsion in the outer Disc, host- 
ing essentially old stars, at odds with, for instance, usual 
spiral arms. 

In our view, t he evidence collected h e re, in Pap- 
I, Pap-II and by iMartinez-Dekado et"al] (l2005al) . to- 
gether with the kinematics of Canis Major stars as de- 
scribed in Pap-III, Pap-I V and the po ssible links with 
the Monoceros Ring (Pap-I , Pap- IV, |P^naimhiaet_alJ 
[ 20051: iMartinez-Delgado et al] l2005bl: iDinescu et al] 120051: 
lcmmet^uTl27)fI t )bl) are more naturally explained by the hy- 
pothesis that Canis Major is the remnant of the disrupting 
dwarf galaxy in a nearly-circular and nearly-planar orbit 
about the centre of the Galaxy. Accepting this hypothesis 
as the one that best fits all the available data we have derived 
several physical characteristics of CMa as a stellar system, 
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Table 2. Fundamental parameters of the main body of the Canis 
Major system. 



{D sun ) 

lo a 
bo b 

(Rgg) 



y 

z 

My 

/V,240° 

A*V,244° 

HWHMi.o., 



7.2 ±1.0 kpc 
244.0° ± 2.0° 
- 6.0° ± 1.0° 
13.1 ± 1.0 kpc 
3.7 ±0.6 kpc 
- 6.4 ± 0.6 kpc 

~ -0.8 kpc 
- 14.4 ± 0.9 mag 
24.0 ± 0.6 mag/arcsec 2 
23.8 ± 0.6 mag/arcsec 2 
~ 2.0 kpc 



a As derived from Fig. 8. Adopting the distance selected sample of 
Fig. 10, l = 241.7° ± 1.0° is obtained (see Sect. 3.2.1). 
" As derived from Fig. 8. 

that are summarised in Table 2. In particular we draw the 
attention of the reader on the following specific results of 
the present analysis: 

• RC stars provide a quantitative description of CMa that 
is in full agreement with those obtained from other inde- 
pendent trac ers (M giants and MS stars, see Pap-I, Pap-II, 
Pap-III and iMartinez-Delgado et alj |2005a'l. In particular 
we obtained a (distance-independent) estimate of the inte- 
grated V magnitude and a (distance-dependent) estimate of 
the central s urface brightness in good agreem ent with those 
obtained bv lMartinez-Delgado et all (l2005aT) from deep op- 
tical photometry of MS stars. We also obtained a distance 
profile (along the line of sight) of CMa that is essentially 
indistinguishable from that obtained in Pap-III from M gi- 
ants. The high degree of self-consistency achieved among the 
various analyses — using different tracers — may be consid- 
ered also as a validation of our criterion of selection of RC 
stars. We confirm, with independent tracers and an indepen- 
dent distance scale with respect to Pap-I and Pap-III, that 
the mean distance of CMa is Dq = 7.2 ± 1.0 kpc, around 
I = 240° . 

• We have detected a clear relation between the mean 
distance to CMa and Galactic longitude within the main 
body. The derived spatial orientation of the system is in good 
agreement with the predictions of the N-body simulation 
presented in Pap-IV, that models CMa as a dwarf galaxy 
being accreted in a planar orbit onto the disc of the Milky 
Way. 

• Using the same dataset (2MASS), the same tracer (RC 
stars) and the same distance scale as in the analysis sum- 
marised above, we showed that the old open clusters AM 2, 
To 2 and Haf 11 are located (in space) within the main body 
of Canis Major. 

In summary, the present study fully (and indepen- 
dently) confirms th e results presented in previous papers 
of this series and in lMartinez-Delgado et al ] (l2005al) : Canis 
Major has the size, the luminosity and the kinematics typi- 
cal of a large dwarf galaxy, seen during the last stages of its 
disru ption within the tidal field of the Milky Way (see also 
ISbordone et al]l2005lh While a Galactic origin for CMa can- 
not be definitely ruled out with the present data, the dwarf 



galaxy hypothesis appears as best suited to fit the overall 
observational scenario. 

While further details emerge — as, for example, the 
distance-longitude gradient — we begin to obtain a clearer 
view of this challenging system. Future studies should try 
to overcome the several observational challenges posed by 
this nearby, low latitude object, to unveil the details of its 
stellar content (age, chemical composition) and to clarify its 
possible connection with the Monoceros Ring. 
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APPENDIX A: CONTAMINATION BY 
DWARFS 

To study the effect of the contamination by dwarfs in our 
color-selected sample of candidate RC stars o n subtracted 
densit y maps, we recur to the Galactic model bv lRobin et alJ 
(2003). This model includes the effect of the Galactic Warp 
and of the disc flaring, and takes into account the off-plane 
position of the Sun. The main advantages of using the R03 
model for our purpose can be summarized as follows: 

• we can obtain (K,J-K) CMDs of the synthetic samples, 
therefore we can apply exactly the same selections and anal- 
ysis adopted for the observed sample; 

• synthetic stars from the R03 models are individually 
classified into different species with the flag cl, hence the 
discrimination of dwarfs from giants is straightforward. In 
particular stars having cl = 3 are giants, those with cl = 4 
are subgiants and those with cl = 5 are dwarfs. Since the 
synthetic CMDs reveal that both cl — 4 and cl — 5 stars 
can fall into our selection color-window, from now on we 
will consider all the stars of class 4 or 5 as dwarfs. 

The obvious disadvantage is that the result will not re- 
flect the real Galaxy but just a model of it. However, in the 
present context, the unavoidable approximations inherent to 
a Galactic model should not represent a major concern. In 
the following we will mainly focus on the magnitude range 
12.0 ^ K 14.0, that is the range enclosing essentially all 
the stars that contribute to the main structures we discuss 
in the present paper, and CMa in particular. 

To allow the possibility of a direct comparison with ob- 
served CMDs (a very useful sanity check) and to limit the 
dimensions of the synthetic catalogues we choose to work 
on 30 deg 2 windows in the latitude range 6° ^ |6| ^ 11°, 
i.e. exactly the same as the observed sample whose CMD is 
shown in Fig. 1. To study the dwarf contamination toward 
different directions of interest, we selected three windows 
at / = 180° ± 3°, I = 241° ± 3° (i.e., the same as Fig. 1), 
and I = 280° ± 3°. We extracted from the model the cor- 
responding samples in the Southern and Northern Galactic 
hemispheres, at the above longitudes, for a total of six cat- 
alogues. The extraction was performed using the large field 
option, that takes into account the variations of the popu- 
lation across the considered fields 4 . To reproduce as well as 
possible the real case we let the model include the effects of 
interstellar extinction according to the default assumptions 
included in the web form. 

From the observed dataset of Fig. 1 we fitted the follow- 
ing relations for the photometric errors in K(J) as a function 
of - not extinction-corrected - K(J) magnitude: 

eK = 0.0236 - 0.00047A: + 3.480 x 10~ 8 exp(K) mag (Al) 

eJ = 0.0169 + 0.00024J + 5.840 x 10~ 9 exp( J) mag . (A2) 

To simulate realistically the effect of photometric errors 
on synthetic stars we computed eK and eJ for each of them 
with the above equations, we randomly extracted a Gaussian 
deviate (8, positive or negative, in a units) for each of them, 
and then we added the resulting synthetic errors 8k X eK 

4 See the web form of the R03 model and the explanations pro- 
vided therein, |http:/ /bison. obs-besancon.fr/modele/| 
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Figure Al. Example of synthetic CMDs extracted from the R03 
model. The symbols and the scale are the same as in Fig. 1. The 
continuous line encloses the color-window 0.45 J — K sC 0.70. 



and 8 j x eJ to their K and J magnitudes, respectively. After 
that, we corrected the synthetic stars for extinction, using 
the Av value associated by the model to each star. In the 
following J and K stand for extinction-corrected magnitudes, 
as in the rest of the paper. 

The above procedure should provide synthetic samples 
as similar as possible to the observed ones, including the ef- 
fect of photometric errors. The comparison of the final syn- 
thetic CMDs for the I = 241° window (Southern and North- 
ern hemisphere) shown in Fig. Al with the observed CMD 
of Fig. 1 demonstrates that this objective has been achieved: 
the overall similarity between observed and synthetic CMDs 
is striking. In Fig. Al we report also the adopted color win- 
dow we use to select RC stars in the observed sample. A 
shift of 0.04 mag in J-K has been applied to all the synthetic 
samples to center the vertical plume of RC stars within the 
selection window, as it is in the observed case. This small 
mismatch in color is likely due to a non-perfect match be- 
tween the theoretical star tracks used by the R03 model 
and real stars (quite a usual occurrence) as well as to tiny 
differences in the adopted extinction laws. In any case the 
following results are essentially unaffected by the small color 
shift applied. In the following we will consider only syntethic 
stars falling in our color-selection window. Hence, when we 
speak of dwarfs we mean stars of class 4 or 5 (R03) having 
0.45 ^ J — K ^ 0.70, and for giants we mean stars of class 
3 having 0.45 ^ J — K ^ 0.70, that is, the stars that we 
would have selected as RC stars in our analysis. 

Before proceeding with our set of tests, Fig. Al deserves 
a few comments. It is quite clear that there is a strong 
asymmetry between the Southern and Northern samples: 
the Southern sample contains many more stars and the dif- 
ference is particularly evident in the RC plume. This is the 
effect of the Warp included in the model. However the com- 
parison of the LFs of the observed (Fig. 1, left panel) and of 
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Figure A2. Fraction of dwarf stars in the total number of color- 
selected RC stars as a function of K magnitude. The three panels 
show Fp for three different directions, in the Northern (open 
squares) and Southern (x symbols) hemispheres. 



the synthetic (Fig. Al, left panel) Southern samples shown 
in the right panel of Fig. 1 demonstrates that (a) the R03 
model provides an excellent overall representation of the ob- 
served distribution of stars in the considered color window 
(except for the total normalization, which appears too large 
by a factor ~ 3), and (b) it still fails to reproduce the strong 
bump in the observed LF around K ~ 13.0, i.e. the putative 
RC of CMa. This comparison is fully consistent with the 
coesistence of the Galactic Warp and an additional stellar 
system, CMa. Finally, the excellent agreement between the 
syntetic and observed LFs outside the K = 13.0 ±0.5 range, 
seems to confirm that our sam ples are essentially complete 
down to K ~ 14.5, as stated in lCutri et al] <l2003l) . 

In Fig. A2 the fraction of dwarfs falling into our color- 
selection as a function of K magnitude (Fd) is shown for the 
Northern (open squares) and Southern (crosses) samples for 
the three selected directions. There are several interesting 
conclusions that can be drawn from this plot: 

• First of all, since the line of nodes of the warp included 
in the R03 model lies exactly toward the I = 180° direction, 
any South-North difference in Fp and/or in the detected 
South-North asymmetries in the I = 180° sample must be 
due to the effect of the ~ 15 pc displacement of the Sun to 
the North of the Galactic Plane. The upper panel of Fig. A2 
shows that difference in Fd is very small, less than 5 per cent 
in any bin for K 12.0 and ~ 1 per cent in the whole 12.0 ^ 
K ^ 14.0 range. Moreover, the South-North difference in 
giants at / = 180° is 5 per cent of that at I = 241° and 2 per 
cent of that at I = 280°. Finally the South-North difference 
in the total number of color selected RC stars at I = 180° 
is 12 per cent of that at I = 241° and 6 per cent of that 
at / = 280°. The conclusion is that the off-plane position 
of the Sun has very little effect both on the contamination 



by dwarfs and on the amplitude of observed South-North 
asymmetries, for the applications considered here. Hence our 
decision to neglect it in the analysis described in this paper 
appears to be fully justified. 

• In the range most relevant for the present study (12.0 ^ 
K ^ 14.0) Fd is quite large, reaching the 80 per cent level 
at K = 14.0, in some cases. 

• The overall degree of contamination is quite similar in- 
dependent of the considered directions, from the Anticenter 
to / = 280°. 

• There are significant differences in the Northern and 
Southern F D toward I = 241° and I = 280°, up to 20 - 30 
per cent in some bins. However Fd is always higher in the 
Northern hemisphere than in the Southern one. It appears 
that, according to the R03 model, N-S differences in the 
contamination by dwarfs are working against the detection 
of real overdensities in the Southern hemisphere. 

The latter result may appear puzzling at a first glance, 
but it can be easily explained by looking into the synthetic 
samples. It turns out that it is due to the presence of a 
significant South-North asymmetry at large Galactocentric 
distances (i.e. the Warp) that is much larger in giants than in 
dwarfs, in the magnitude ranges of interest. For example, in 
the range 12.0 < K < 14.0, at I = 241° the South-North ra- 
tio of dwarf density is p D {S)/p D (N) = 29096/24371 = 1.19 
while for giants p G (S)/p G (N) = 28684/7572 = 3.79. There- 
fore, even if there are more dwarfs in the South than in 
the North, the imbalance in giants is much larger and 
the final Fd is consequently larger in the North (Fd =76 
per cent) than in the South (Fd =50 per cent). At 
/ = 280°, p D (S)/p D (N) = 35427/31333 = 1.13, while 
Pg(S)/pg(N) = 70305/22067 = 3.18, in the same range 
of magnitudes. This implies that in presence of a distant 
South-North asymmetry (as the Warp, as in the R03 model 
we are presently considering, or CMa, or both), the observed 
imbalance in dwarfs is small in this magnitude range, such 
that the contribution by dwarfs essentially disappears when 
Southern and Northern counts are subtracted. 

The above results drive us to the two key questions we 
have to answer to really check the reliability of the results 
of the present paper: 

(i) What is the contribution of dwarf contamination to 
subtracted densities! 

(ii) Is our observable, the subtracted density of color se- 
lected RC stars p(S) — p(N) , a good tracer of the real over- 
density of giants? In other words, is it a good estimator of 
the quantity we are interested in (pc(S) ~ pc(N))? 

Fig. A3 answers to the first question. There we plot the 
fraction of the total S-N density difference that is due to 
dwarfs 

Ap D = p D (S) - pp(N) 
Ap T ~ p(S) - p(N) 

as a function of magnitude, for the three considered di- 
rections. The subscript T stands for Total=Dwarfs+Giants. 
The ratio of the total Southern-to-Northern density is also 
reported within each panel. This plot shows that in absence 
of a real South-North asymmetry (I = 180°, p(S)/rho(N) ~ 
1) the contribution is quite significant (> 40per cent) for 
K > 12.0. On the other hand, toward directions that cross 
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Figure A3. Contribution of the South-North imbalance in dwarf 
stars to the total South-North density difference as a function of 
magnitude, for three different directions. p(S)/p(N) is the ratio 
of color-selected RC stars in the range 12.0 ^ K ^ 14.0, and 
roughly describes the degree of total South-North imbalance in 
each give direction. 



Figure A4. Correlation between the total subtracted density of 
color-selected RC stars and the real subtracted density of giants. 
The continuous line passes through the origin and has slope 1.0, 
while the dotted line is the best linear fit to the data (see the 
equation in the upper left corner of the plot). The directions cor- 
responding to the three plotted points are indicated in the upper 
x axis. 



the Warp asymmetry (I = 241°, 280°, p(S)/rho(N) > 1.8) 
the dwarf contribution remains lower than 30 per cent every- 
where, for K ^ 14.0. Hence, even in a regime of strong con- 
tamination by dwarfs, the effects on South-North subtrac- 
tion are quite limited for K ^ 14.0 and in the presence of an 
overdensity like the Warp and/or CMa. On the other hand, 
our subtracted maps may be significantly affected (or even 
dominated, in some cases) by dwarf contamination around 
/ = 180° and for K > 14.0. It is very interesting to note that 
in those regions we see a low-significance and smooth South- 
North residual (Fig. 3) that completely disappears when we 
subtract by the rescaled Northern density map (Fig. 5) . On 
the other hand all the structured overdensities we deal with 
in this paper lie in the region where dwarfs provide a mi- 
nor contribution to the subtracted densities (I > 200° and 
K ^ 14.0). Within this context, there is another interesting 
point worth considering. Since the contribution of dwarfs 
to our subtracted densities grows continuously with magni- 
tude even for K ^ 14.0, if our maps were seriously affected 
by dwarf contamination we would expect to see a corre- 
sponding growth of density with heliocentric distance. On 
the contrary, in the direction of CMa we see both a growth 
and a decrease of the density (see Figs. 3 and 5). The dis- 
tance profile of Fig. 11 shows that toward I ~ 240° the 
subtracted density reaches a maximum at Dq — 7.2 kpc, 
corresponding to K = 12.8, and then essentially falls to zero 
at Dq ~ 10.0 kpc, corresponding to K — 13.5, a completely 
different profile with respect to what is expected from a dis- 
tribution of Galactic dwarf stars. This occurrence strongly 
suggests that dwarfs cannot have a major role in determining 
the observed distribution of the CMa overdensity. The same 
argument is valid also when considering the trend of con- 



tamination of subtracted densities as a function of Galactic 
longitude. Since the degree of contamination grows toward 
I — 180°, if dwarfs were dominant contributors to the sub- 
tracted density we should see the maximum density toward 
/ = 180°. On the contrary, the density in that direction is 
low and smooth, while it is much higher and more strongly 
peaked toward I ~ 240° . 

Fig. A4 tries to answer question (ii), above. This is sim- 
ply a plot of the S-N density difference in the total number of 
stars selected in color as RC candidates with our technique 
(i.e., our observable Apr), versus the actual difference in the 
S-N density of giants (i.e. the quantity we want to estimate, 
Apa), in the range 12.0 ^ K ^ 14.0, for the three directions 
we are considering. The correlation between Apr and ApG is 
very strong and essentially linear. Moreover the cofneient of 
the linear best-fit model is near unity, hence our observable 
appears as an excellent estimator of the true South-North 
density difference in giants, at least in the presence of an 
asymmetry like the warp embedded in the R03 model. 

In summary, based on the above experiments with the 
R03 Galactic model, we can conclude that: 

• The neglection of the displacement of the Sun with re- 
spect to the Galactic plane has very little effect on our anal- 
ysis, and can be safely ignored. 

• While our samples of color-selected RC stars may be 
heavily contaminated by dwarfs, the effects of the contami- 
nation are much lower on the subtracted-density maps. The 
contribution of dwarfs to the "observed" Apr is lower than 
30 per cent for K ^ 14.0 toward I = 240° and I = 280°, and 
much lower than this for K ^ 13.0. Note that these num- 
bers take into account only known Galactic components, the 
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presence of an extra asymmetry in the considered range of 
distances (i.e. CMa, see Fig. 1) would further reduce the role 
of dwarfs. 

• The observable we adopt (Apr) appears to trace very 
well the true South-North imbalance of giants under the 
realistic conditions provided by the R03 model. 

It is conceivable that the CMa overdensity we observe 
using color-selected putative RC stars may be in fact pro- 
duced by an unexpected asymmetry in dwarfs. However 
the dwarfs contributing to our counts in the range 12.0 ^ 
K ^ 14.0 have heliocentric distances 0.0 kpc ^ Dq ^ 
4.5 kpc. Their distance distribution shows a strong peak at 
Dq ~ 1.0, the mean is (-Dq) = 1.0 kpc and the standard de- 
viation is a = 0.7 kpc, and 89 per cent have Dq < 2.0 kpc. 
Hence a strong and very compact clump of dwarfs 5 would 
have to be present in the Southern hemisphere, at ~ 1.0 kpc 
from the Sun toward I — 240°, a feature that is obviously 
not present in the R03 model. As far as we know there is no 
report in the literature of such a local overdensity. Finally, 
it should be recalled that the position and size of the RC 
overdensity studied here matches very well those derived for 
CMa with other tracers (M giants and MS stars) that are 
much less affected by contamination from spurious sources. 



5 Actually it must be significantly more compact than the dis- 
tribution of RC stars we attribute to CMa, since the selected 
dwarfs would span a much larger range of absolute magnitudes, 
with respect to RC stars. 



